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I.  OVERVIEW 

A.  Introduction 

The  ground  movement  routines  documented  in  this  report  are  designed 
to  be  part  of  the  STAR  (Simulation  of  Tactical  Alternative  Responses)  ground- 
air  combat  model  developed  at  the  Naval  Postgradute  School.  STAR  is  a  high 
resolution  computer  simulation  model  written  in  the  SIMCRIPT  II. 5  simulation 
language.  The  ground  movement  routines  are  responsible  for  updating  the  po- 
sitions of  vehicles,  (both  combat  and  support)  which  are  moving  on  the  ground 
during  the  simulation.  Separate  routines  will  be  used  for  air  movement. 

The  movement  routines  are  designed  to  allow  a  variety  of  different 
modes  of  movement  and  to  interface  with  terrain  representations  of  varying 
degrees  of  detail.  Thus  if  movement  is  an  important  facet  of  a  study,  the 
terrain  and  vehicle  mobility  data  can  be  incorporated  at  a  high  level  of 
detail,  while  if  mobility  is  felt  to  be  of  subsidiary  importance,  a  less 
detailed  representation  can  be  used. 

The  movement  routines  in  the  report  were  first  developed  for  use  in 
the  updated  STAR  Battalion  model  (ca.  spring  1979),  and  have  been  incorporated 
into  the  expanded  STAR  Brigade  model  which  is  currently  under  development. 

The  remainder  of  the  first  section  of  the  report  will  review  the  basic 
assumptions  of  the  movement  model  and  will  display  the  various  modes  in  which 
it  can  be  used.  Section  II  will  detail  the  data  arrays  which  the  movement 
model  requires.  Section  III  will  discuss  the  ground  movement  subroutines  in 
detail,  and  Section  IV  will  consider  the  interfacing  of  the  movement  model  with 
the  rest  of  STAR.  Finally,  Section  V  will  briefly  indicate  some  aspects  of 
movement  which  are  not  developed  in  this  report.  Discussion  of  other  parts  of 
the  STAR  model  can  be  found  in  references  [1]  thru  [4]. 


B.  Basic  Assumptions  of  the  Movement  Model 

The  ground  movement  model  of  STAR  moves  vehicles  between  posi  tions 
along  routes  which  are  made  up  of  straight  line  segments.  At  any  time  in  the 
simulation  the  STAR  model  can  request  a  location  update  for  any  vehicle. 
Movement  in  the  model  is  by  individual  vehicle.  Movement  control—the  decision 
of  whether  to  move  and,  if  so,  where--is,  however,  usually  by  some  higher  level 
unit.  Thus  the  relationship  of  routes  to  positions  is  organized  by  platoons. 

The  origins  and  destinations  within  the  model  are  represented  by 
Position  areas.  A  position  area  is  a  collection  of  individual  element  X,  Y 
coordinates  (e.g.  a  company  team  defensive  position,  or  the  starting  position 
for  an  attacking  battalion).  The  simulation  model  will  support  an  arbitrary 
number  of  position  areas. 

Within  each  position  area,  the  individual  element  positions  are  or- 
ganized by  platoon.  Movement  from  one  position  area  to  another  is  along 
routes  which  are  assigned  for  each  platoon.  Thus  a  fundamental  assumption 
embedded  in  the  data  structure  for  the  movement  model  is: 

If  two  vehicles,  both  in  platoon  i,  are  to  move  from 

position  area  A  to  position  area  B,  they  must  use 

the  same  route. 

As  we  shall  see  shortly,  this  does  not  require  that  they  follow  ex- 
actly the  same  path  (because  of  the  formation  offsets),  but  essentially  they 
will  move  in  the  same  general  location.  Also,  since  movement  is  by  individual 
vehicle,  we  are  not  required  to  always  move  the  entire  platoon  at  the  same 
time  along  the  same  route,  but  often  this  will  be  the  case. 
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Movement  formations  may  be  used  to  ensure  that  a  platoon  moving 
along  a  route  will  do  so  in  an  organized  fashion  (e.g.  to  achieve  a  line  for- 
mation in  the  final  assault  phase  of  an  attack:  formation  offsets  displace 
the  platoon  members  left  and  right  from  the  route  they  are  following). 
Movement  for  each  vehicle  is  continuous  in  the  sense  that: 
i)  Any  x,  y  coordinate  on  the  battlefield  is  a  possible 

location  for  a  unit, 
ii)  Arbitrarily  small  (or  large)  movement  increments  can  be 
requested  (e.g.  move  a  unit  for  2.2  seconds), 
iii)  The  movement  speed  for  a  unit  is  continuous  —  if  a  unit 
wants  to  change  speed,  it  can  do  so  only  gradually,  lim- 
ited by  bounds  on  possible  acceleration  and  deceleration, 
iv)  Acceleration  however,  is  not  continuous.  Typically  if 
a  unit  wishes  to  accelerate  it  will  do  so  at  the  maximum 
allowed  acceleration  rate  until  the  new  speed  is  reached, 
and  then  will  move  at  that  new  (constant)  speed  until  it 
is  time  to  change  speed  again. 
Limitations  on  speed  and  acceleration  may  be  derived  for  each  vehicle  from 
actual  terrain  data  (for  high  resolution  movement)  or  may  be  set  to  reason- 
able values  for  each  vehicle  throughout  the  simulation  (for  lower  resolution 
movement) . 

C.  Capabilities  of  the  Model 

This  section  summarizes  the  various  modes  in  which  the  movement  model 
may  be  used.  These  modes  relate  to  the  choice  of  origin,  destination,  route, 
and  formation  for  the  movement. 
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1.  MOVEMENT  IN  PREDETERMINED  DIRECTION 

For  simplistic  analyses,  the  model  can  be  used  to  analyze  an  attack 
where  each  element  is  given  a  starting  position  and  a  direction  of  movement. 
No  route  need  be  specified  as  movement  is  always  along  the  specified  direction. 

2.  MOVEMENT  TO  SPECIFIED  POSITION  IN  A  POSITION  AREA 

Given  the  current  location  of  a  vehicle,  the  move  model  can  move  it 
directly  along  a  straight  line  (no  route)  to  a  specified  position  within  a 
given  position  area.  This  mode  would,  for  example,  allow  direct  movement 
from  one  position  area  to  another. 

3.  MOVEMENT  TO  BEST  POSITION  IN  A  POSITION  AREA 

Mode  3  is  similar  to  mode  2,  except  that  the  model  scans  all  possible 
positions  in  the  vehicle's  platoon  area  (in  the  given  destination  position  area) 
and  selects  the  best  position  which  is  not  currently  occupied. 

4.  MOVEMENT  TO  DESIGNATED  ROUTE 

Given  any  current  vehicle  position  and  a  choice  of  route,  the  move  model, 
will  move  the  vehicle  from  its  position  onto  the  route. 

5.  MOVEMENT  ALONG  DESIGNATED  ROUTE 

Given  a  position  on  a  route,  the  model  can  continue  along  the  route — 
(this  is  probably  the  most  frequently  used  mode). 

6.  ROUTE  SELECT 

Given  origin  and  destination  position  area  numbers  the  model  will  select 
the  route  to  be  used  and  then  proceed  as  in  modes  4  and  5. 

7.  MOVEMENT  WITH   FORMATION  OFFSET 

In  modes  5  and  6,  the  movement  path  along  the  route  may  be  offset  from  the 
route  to  put  different  platoon  members  in  different  relative  positions  within 
a  movement  formation.  The  current  implementation  has  platoon  formations  only. 
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8.  FORMATION  SPECIFIED  BY  TERRAIN 

Normally  a  movement  formation  will  be  tactically  determined.  Sometimes, 
however,  the  terrain  forces  a  particular  formation  (e.g.  column  for  crossing 
a  bridge).  The  move  model  can  store  and  automatically  implement  formation 
changes  which  are  required  by  the  terrain. 

9.  STOP  ALONG  ROUTE 

Upon  command,  the  move  model  will  stop  a  vehicle  which  is  moving.  This 

is  useful,  for  example,  if  an  attacking  force  decides  to  enter  a  hasty  defense 

because  of  attrition  levels. 

10.  REVERSE 

Upon  command,  the  move  model  will  reverse  the  previous  direction  of 

movement  of  a  vehicle  and  cause  it  to  return  to  the  position  from  which  it  began 

This  feature  could  be  used  if  a  thwarted  attacking  force  has  to  retreat. 

These  modes  are   automatically  combined  in  certain  circumstances.  For 

example,  a  command  to  move  from  area  A  to  area  B  might  automatically  invoke 

modes 

i     6     (to  select  the  route) 

4     (to  get  onto  that  route  from  area  A) 

5  and  7  (to  move  along  the  route  in  formation) 
(maybe  8) 

3    (to  get  into  the  best  available  position 
upon  arrival  in  area  B) 

The  control  commands  which  determine  the  modes  used  will  be  discussed  later 

in  the  report. 
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II.   DATA  REQUIREMENTS 

In  this  section  the  data  required  to  support  the  ground  movement  model 
is  specified.  As  indicated  in  the  previous  section  there  are  several  different 
modes  of  movement  which  can  be  used.  Depending  on  the  mode  used,  the  data 
requirements  may  vary.  The  general  model  which  allows  movement  between  arbi- 
trary position  areas  along  routes  and  using  formation  offsets  requires  that  all 
of  the  following  be  available. 

A.  POSITION 

POSITION  is  a  3-dimensional  ragged  array  which  gives  the  X  and  Y 
coordinates  of  potential  vehicle  locations  (e.g.  prepared  defensive  positions, 
starting  attack  positions).  Also  included  is  an  orientation  angle,  0  ,  for 
each  position  to  allow  vehicles  in  the  defense  to  select  their  sector  of  re- 
sponsibility. 9  is  measured  in  radians  counterclockwise  from  East.  Positions 
are  organized  by  platoons,  and,  for  each  platoon,  by  the  position  areas  which 
that  platoon  might  occupy. 

The  array  POSITION  (I,J,K)  has  subscripts 

I  =  platoon  number  (I  =  1,...,PNUM) 

J  =  counter  for  areas;  J  =  1,...,  number  of  areas  this  platoon 

might  occupy 

(J  may  be  different  for  different  platoons,  I.) 

For  the  Jth  position  area  for  Platoon  I,  the  X,  Y,  0  values 

are  organized  as  follows.  Typically  all  areas  for  a  given  platoon  will  have 

as  many  vehicle  positions  as  there  are  elements  in  the  platoon. 
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K  =  1     Position  Area  ID  number 

2  X 

3  Y  \        x>  y»  Q   f°r  first  vehicle  in  this 

1  /  platoon  in  this  position 

4  0]  ■  area 

5  X2 

6  Y?  \  second  position  for  this 

platoon  in  the  position 

02  area 


3n-l 

X 
n 

3n 

Y 
n 

3n+l 

en 

last  position  for  this 
platoon  in  the  position 
area 


If  the  movement  model  is  to  be  used  in  mode  3--select  best  position 
in  platoon  area-- the  individual  vehicle  positions  1,  ...  ,  n  are  assumed  to 
be  arranged  in  priority  order.  (X,Y,0,  is  the  best  position,  X?Yp02  is  next 
best,  etc.)  This  mode  can  be  important  in  moving  to  subsequent  defensive 
positions  in  the  active  defense.  If  only  part  of  a  platoon  survives  to  reach 
the  new  position  area,  we  want  them  to  occupy  the  vehicle  positions  which  have 
the  best  defensive  characteristics. 

The  POSITION  array  is  initialized  by  the  following  segment  of 
Simscript  code  in  the  RES. MOVE  program  given  in  Figure  1. 

To  illustrate  the  POSITION  array  and  other  arrays  to  be  defined  in  this  section, 
consider  the  following  condensed  and  simplified  movement  scenario. 

Platoon  1  has  3  vehicles  and  will  fight  in  place  in  position 
area  37--no  movement  required. 

Platoon  2  has  2  vehicles.  It  begins  in  position  area  37  and  may 

move  either  to  position  area  16  or  to  position  area  42  along  designated  routes. 
(See  Figure  2.)  -  7  - 


1  ROUTINE  RES. MOVE 

2  DEFINE  I,  J.  K,  N,  ID.  NE.  Nfl.  NM  fiS  INTEGER  VflRIflBLES 

3  USE  UNIT  5  FOR  INPUT 

4  USE  UNIT  6  FOR  OUTPUT 

5  LET  MRX.DIST.INCR  =  50. 

6  LET  FOR.CHG.INT  =  200. 

7  RESERVE  POSITION  (x, *, *)  AS  PNUM  BY  k   '  "PNUM  IS  NUMBER  OF  PLATOONS 

8  FOR  I  =  1  TO  PNUM 

9  DO 

10  READ  ID.  NE.  Nfl   '  "PLT  NO.,  NO.  ELEMS  IN  PLT.  NO.  AREAS  USED  BT  PLT 

11  IF  ID  NOT  EQUAL  TO  I  PRINT  1  LINE  WITH  I  AS  FOLLOWS 

12  XXXXX  INPUT  DATA  SEQUENCE  ERROR  IN  POSITION  ARRAY  FOR  PLT  *****  XXXXX 

13  ALWAYS 

14  RESERVE  POSI T ION (I , *. *)  AS  NA  BY  3*NE  +  1 

15  FOR  J  =  1  TO  NA 

16  DO 

17  FOR  K  =  1  TO  3*NE  +  1 

18  READ  POSITIGNtl.  J.K)   "ALL  POSITION  DATA  FOR  THIS  PLATOON 

19  FOR  K  =  14  TO  3*NE+1  BY  3 

20  LET  POSITION (1. J. Kl  =  POSITION  (I .  J,  K)  /RADIAN.  C   "DEG  TO  RADIANS 

21  LOOP 

22  LOOP    "END  OF  INITIALIZATION  FOR  POSITION  ARRAY 


Figure  1,  Initialization  of  POSITION  Array 
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Position  Area  37 


PLT  1 


o  ^ 


Route  2 


Figure  2.  Simplified  Movement  Scenario 
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Input  data  for  the  POSITION  array  consistent  with  this  movement 

scenario  follows: 

1  3    1  (platoon  #1  has  3  vehicles  and 

may  occupy  1  position  area) 


37  X1  Y]  e1  X2  Y2  e2  X3  Y3  93      ^the  single  area  is  area  37  and 

within  area  37  the  3  vehicles  have 
X,  Y,  and  9  coordinates) 


(platoon  #2  has  2  vehicles  and 
may  occupy  3  position  areas) 


37  X4  Y4  e4  X5  Y5  95  The  3  areas  are  areas  37,  16  and 

16  Xc  Yc  dc   X,  Y,  e_  42.  Each  has  2  vehicle  positions 

6  6  6  7  7  7 

42  X8  Yg  98  Xg  Yg  9g 


B.  MOVE. DATA 

The  MOVE. DATA  array  is  used  to  select  a  route  between  two  position 
areas  for  a  given  platoon.  The  data  is  organized  by  platoon.  For  each  pla- 
toon, MOVE. DATA  contains  a  list  of  triples 

where  A,,  A?  are  area  numbers  and  R  is  the  number  of  the  route  connecting 
A,  to  A?  .  For  each  triple  it  is  assumed  that  A,  is  less  than  A?,  and 
that  route  R  starts  at  A,  and  ends  at  A«  .  If  a  vehicle  wants  to  move 
from  area  A-  to  area  A,,  then  it  will  traverse  route  R  backwards. 
The  MOVE. DATA  array  has  subscripts 

I  =  platoon  number  (I  =  1,...,  PNUM) 

J  =  1,...,  3  *  number  of  area  pairs  for  this  platoon. 
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For  platoon     I,   the  data  is  organized  as   follows 
J  =  1         A- 


1 


Area  Numbers 


2    A, 


3  RT(1  ->  2)   Route  from  Al  to  A2  for  platoon  I 

4  A. 


5    A, 


} 


Area  Numbers 


6    RT(3  -*■  4)   Route  from  A3  to  A4  for  platoon  I 


The  (A, ,A„)  pairs  in  the  list  are  assumed  ordered  lexicographically  in 
increasing  order  to  aid  the  search  process  in  route  selection  (e.g.  if  the 
area  pairs  to  be  connected  are 

(1,2),     (1,4),     (2,4),     (3,2) 

then  their  order  in  the  list  is  as  written  above). 

Note  that  different  platoons  may  use  different  routes  between  the 
same  area  pairs.  Also,  that  only  those  area  pairs  which  are  used  by  a  platoon 
will  be  included  in  its  section  of  the  MOVE. DATA  Array. 

MOVE. DATA  is  initialized  in  the  RES. MOVE  program  by  the  code  segment 
in  Figure  3. 

A  typical  set  of  input  data,  consistent  with  the  previous  POSITION 
data  follows:   (again  PNUM  =  2  =  #  of  platoons) 


1  0 

2  2 

16   37   2   37   42   1 


platoon  1  has  only  one  area,  so  it 
will  never  move—zero  routes  used, 
platoon  2  has  2  movement  possibilities 


Platoon  2  will   use  route  2  to  get  from  area  16  to  area  37  or  route  1   to  get 
from  37  to'  42 . 
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23  RESERVE  MOVE.  DATA  (x,  x)  AS  PNUM  BY  x   "PNUM  IS  NUMBER  OF  PLATOONS 

2U  FOR  I  »  1  TO  PNUM 

25  00 

26  READ  ID,  NM     "PLATOON  NO.,  NO.  OF  ROUTES  USED  BY  THIS  PLT 

27  IF  ID  NOT  EQUAL  TO  I  PRINT  1  LINE  WITH  I  AS  FOLLOWS 

28  XXXXX  INPUT  DATA  SEQUENCE  ERROR  IN  MOVE. DATA  ARRAY  FOR  PLT  xxxxx  XXXXX 

29  ALWAYS 

30  IF  NM  EQUALS  0    CYCLE 

31  ELSE 

32  RESERVE  MOVE. DATA  (I . x)  AS  3xNM 

33  FOR  J  =  1  TO  3xNM 

34  READ  MOVE.DATAtl. J) 

35  LOOP    "END  OF  INITIALIZATION  FOR  MOVE. DATA  ARRAY 


Figure  3.     Initialization  of  MOVE. DATA  Array 
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There  could  also  be  a  route  from  16  to  42  but  in  this  example  we 
have  (arbitrarily)   decided  that  the  simulation  scenario  does  not  require  it. 
If  it  were  included,  the  lexicographic  ordering  would  require  that  it  be  be- 
tween the  above  two  triples.     Note  also  that  the  route  runs  from  16  to  37  even 
though  the  platoon  will   be  moving  from  37  to  16. 
C.      ROUTE. DATA 

The  ROUTE. DATA  array  contains  the  description  of  each  of  the  routes 
used  by  the  model.     A  route  is  a  sequence  of     X,  Y     coordinates  called  move- 
ment control   points   (MCPS).     Between  MCPS  the  model   plots  straight  line  route 
segments. 

Subscripts   for  the  array  ROUTE.DATA(I ,  J)  are 

I  =  route  number,     I  =  1,..., number  of  routes 

for  each  route     I,   the     J     subscript  indexes  a  list  of  coordinate  pairs  of 
MCP ' s . 

Along  with  each  MCP  coordinate  pair  is  stored  a  platoon  formation 
code.     If  this  code  is     0     it  means  any  formation  can  be  used  on  this  route 
segment.     Otherwise  the  indicated  formation  must  be  used. 

J  =  1         X. 


Al 

2 

Yl 

3 

form  code. 

4 

h 

5 

Y2 

6 

form  code2 

} 

} 


First  MCP 

Refers  to  segment  between  MCP's  1   and  2 

Second  MCP 

Refers  to  segment  between  MCP's  2  and  3 


The  number  of  MCP's  can  be  different  for  different  routes,  and  many  different 

platoons  may  use  the  same  route.     ROUTE. DATA  is   initialized  in  the  RES. MOVE 

program  by  the  code  segment  in  Figure  4. 
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36  READ  N  "NUMBER  OF  ROUTES  TO  USE 

37  RESERVE  ROUTE. DATA  (*, *)  AS  N         BY  * 

38  FOR  I  =  1  TO  N 

39  DO 

40  REAO  ID.  NM      "ROUTE  NUMBER,  NO.  OF  MVMT  CONTROL  POINTS  ON  THIS  ROUTE 

41  IF  ID  NOT  EQUAL  TO  I  PRINT  1  LINE  WITH  I  AS  FOLLOWS 

42  XXXXX  INPUT  DATA  SEQUENCE  ERROR  IN  ROUTE. DATA  ARRAY  FOR  ROUTE  *x**x  XXXXX 

43  ALWAYS 

44  RESERVE  ROUTE.  DATA  (1 ,  m)  flS  3*NM 

45  FOR  J  =  1  TO  3*NM 

46  READ  ROUTE. DATA  (I. J) 

147  LOOP      "END  OF  INITIALIZATION  FOR  ROUTE. DATA  ARRAY 


Figure  4.      Initialization  of  ROUTE. DATA  Array 
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A  typical  data  input  is  as  follows  (again  consistent  with  our  previous 
scenario) . 

number  of  routes  in  the  array 

route  #1  has  3  MCP's 

(3  MCP  coordinates,  formation  1  must 
be  used  between  MCP2  and  MCP  3) 

route  #2  has  2  MCP's 


2 

1 

3 

Xl 

Yl      ° 

h 

Y2     ' 

X3     Y3     ° 

X,  Y4  0  X5 


0 


(2  MCP  coordinates,  no  formation 
requirement) 

(Note  that  the  X,Y  coordinates  here  are  different  numbers  from  those  in 

the  POSITION  array,  even  though  we  have  used  the  same  notation  for  both) 


D.   FORM.  OFFSET 

The  FORM. OFFSET  array  contains  information  on  the  relative 
position  of  vehicles  in  a  formation  expressed  as  offsets  from  the  route. 
The  offsets  are  assumed  to  be 


MCP 


MCP 


route  segment  and 
direction  of  movement 


Ax   positive  in  direction  of  movement,  and 
Ay   positive  left  of  the  route. 
The  F0RM.0FFSET(I,J)  array  has  subscripts, 


I  =  1, 
J  =  1, 


,  number  of  formations 


,  2  *  number  of  places  in  the  formation  (may  be 
different  for  different  I) 
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and  for  formation    I,  the  offsets  are  organized  as 

J  =  1         Ax,  offsets  for  first  place 

2  Ay,        in  formation 

3  ax?        second  place 


4    Ay, 


etc 


FORM. OFFSET  is  initialized  by  the  final  segment  of  code  in  RES. MOVE  given 
in  Figure  5. 


A  typical 

input  data  set  is 

2 

1         2 

0         0 

-50       0 

I 


number  of  formations 

offsets  for  2  places  in  formation  1; 
place  2  is  50  meters  behind  place  1; 
both  are  on  the  route  in  column. 


0       -150         0       -50       0       50       0       150       (a  4  place  line  formation  with 

100  meter  spacing) 

E.     Vehicle  Attributes 

Each  vehicle  in  the  STAR  simulation  is  a  temporary  entity  which  has 
numerous  attributes  attached  to  it.     The  attributes  which  are  of  interest  to 
the  move  routine  are  listed  below. 

MV. STATE  the  primary  control   variable  for  initiating  and  stopping 

movement,   possible  values  and  their  meanings  are: 

0  -  in  position.     Do  not  move. 

1  -  want  to  start  movement  from  one  position  area  to 

another—do  a  route  selection  and  start  to  move 

2  -  continue  movement  along  a  previously  selected  route 

3  -  stop  along  route  (e.g.   stop  to  fire) 
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48  READ  N  "NUMBER  OF  MOVEMENT  FORMATIONS 

49  RESERVE  FORM. OFFSET (*, *)  AS  N        BY  * 

50  FOR  I  =  1  TO  N 

51  DO 

52  READ  10, NM      "FORMATION  NO.,  NO.  OF  PLACES  IN  THAT  FORMATION 

53  IF  ID  NOT  EQUAL  TO  I  PRINT  1  LINE  WITH  I  AS  FOLLOWS 

54  XXXXX  INPUT  DATA  SEQUENCE  ERROR  IN  FORM. OFFSET  ARRAT  FOR  FORMATION  *****  XXXXX 

55  ALWAYS 

56  RESERVE  FORM. OFFSET  (I , *)  AS  2*NM 

57  FOR  J  =  1  TO  2*NM 

58  READ  FORM. OFFSET  (I.J) 

59  LOOP      "END  OF  INITIALIZATION  FOR  FORM. OFFSET  ARRAY 

60  RETURN 

61  END 


Figure  5.      Initialization  of  FORM. OFFSET  Array 
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AREA. START 


AREA. END 


ROUTE 


NEXT. MCP 


DIR.ON.RT 


X.  CURRENT 


Y. CURRENT 


Z. CURRENT 


SPD 


T.SPD 


DIR.OF.MVMT 


PLT 


4  -  next  position  has  been  reached--stop. 

5  -  final   position  has   been  reached—never  move  again 

The  origin  and  destination  position  area 
numbers  for  the  movement  of  a  vehicle. 

The  route  number  along  which  the  vehicle  moves. 
0  if  not  using  routes. 

MCP  number  on  designated  ROUTE  toward    which  the 

vehicle  is  now  moving. 

0     if  end  of  route  has  been  reached. 

0  if  vehicle  is  moving  in  direction  of  increasing 
MCP  numbers  along  route 

1  if  vehicle  is  traversing  MCP's  in  decreasing  order 
(backwards  along  route) 

location  of  vehicle 

at  most  recent 

movement  update 

speed  of  vehicle  at  end  of  most  recent 
movement  update 

simulation  time  at  which  most  recent  movement 
update  ended.  (Time  SPD  was  last  set.) 

angular  direction  of  movement  (measured  in 
radians  from  East) 


the  platoon  to  which  the  vehicle  belongs 
POS. IN. PLT. AREA  position  number  indicating  which  position  in  the 

POSITION  array  this  vehicle  is  occupying  or  plans  to 
occupy.  (Zero  while  on  route  if  best  position  is  to 
be  chosen  on  arrival  in  position  area.) 

FORM. CODE      formation  number  for  the  platoon  (0  if  not  in  formation; 
then  vehicle  moves  along  route  without  offset). 

FORM. POS       number  indicating  which  place  in  the  formation  this 
vehicle  should  occupy. 
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F.     Set  Membership 

The  vehicle  temporary  entities  may  belong  to  several   sets.     For 
movement  control   purposes  the  most  important  set  is   the  PLT.UNIT  set  which 
is  owned  by  a  PLATOON. LEADER.     Since  position  areas,  routes,  and  formations 
are  organized  by  platoon,  it  is  frequently  necessary  to  reference  a  vehicle's 
platoon  and  the  other  elements  in  that  platoon  by  looping  through  the 
PLT.UNIT  set. 
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III.     MODEL  SUBROUTINES 

In  this  section  we  review  the  subroutines  which  are  used  for  execu- 
ting ground  movement  commands  in  the  STAR  model.     There  are  several   quite 
brief  "utility"   routines,  and  one  quite  lengthy  MOVE  routine.     For  each  rou- 
tine we  list  the  local   variables,   global   variables,   routines  called,  events 
scheduled,  a  brief  description  of  the  code,  and  information  on  how  to  use  the 
routine. 

A.     Routine     INIT.POS 

The  INIT.POS  routine  selects  the  initial   position  of  a  given 
vehicle  from  the  POSITION  array. 
Input  Argument 

VEH 
Local   Variables 


pointer  of  the  vehicle  to  be  positioned 


I,   J,   K 
Global   Variables  Used 

POSITION 
Vehicle  Attributes  Used 

PLT 

POS.IN.PLT.AREA 

AREA. START 

X. CURRENT 
Y. CURRENT 

Z.  CURRENT 

DIR.OF.MVMT 

PRI.DIR 

Routines  called 

BEST. P0S( VEH) 
ELEV 


) 


array  subscripts 

array  of  positions 

platoon  number 

position  number  for  the  vehicle 

the  area  in  which  to  place  this  vehicle 

resulting  X  and 

Y  coordinates  from  POSITION 

elevation  coord  from  ELEV 

both  set  to  resulting  e 

orientation  angle  from  POSITION 
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Events  scheduled 

none 

Code  -  See  Figure  6. 

Brief  Description 

Line  5  IF  POS.IN.PLT.AREA  is  zero,   then  routine 

BEST.POS  is  called  to  select  the  best 

available  position  (this  call   sets 

POS.IN.PLT.AREA  to  a  nonzero  value). 

Line  8  scans  the  list  of  position  areas  for 

this  platoon  until   a  match  is  found  with 
AREA. START 

Lines  9-13  then  select  the  desired  position  from 

this  area  and  set  the  vehicle's     X     and 
Y     coordinates  to  this  position.     The 
vehicle's  movement  and  search  orientations 
are  also  set. 

Line  14  calls  the  ELEV  routine  to  set  the 

vehicle's     Z     coordinate  from  the 
terrain. 

Use 

Enter  with 

i )  vehicle  pointer 

i i )  PLT 

iii)  AREA. START 

iv)  (optional)  POS.IN.PLT.AREA 

On  Exit,  routine  has  set 

i)     POS.IN.PLT.AREA  (if  zero  on  entry) 
ii)      XCURRENT,   Y. CURRENT, Z. CURRENT 
iii)     DIR.OF.MVMT,   PRI.DIR 
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1  ROUTINE  FOR  INIT.POS  CVEH) 

2  "  SELECTS  INITIAL  POSITION  FOR  AN  ELEMENT 

3  DEFINE  VEH.  I,  J,  K  AS  INTEGER  VARIABLES 

4  LET  I  =  PLT  (VEH) 

5  IF  POS. IN. PLT. AREA  (VEH)  EQUALS  0  CALL  BEST. POS (VEH) 

6  ALWAYS 

7  LET  K  =  POS. IN. PLT. AREA (VEH)  *  3 

8  FOR  J  =  1  TO  DIM. F  (POSITION  (I,*.*))  WITH  POSITION  (I . J, 1 )  EQUALS  AREA. START  (VEH) 

9  DO 

10  LET  X. CURRENT (VEH)  =  POSITION (I , J. K-l) 

11  LET  Y. CURRENT  (VEH)  =  POSITION  (I . J. K) 

12  LET  DIR.OF.MVMT  (VEH)  =  POSITION  (I . J, K+l) 

13  LET  PRI.DIR(VEH)  =  DIR.OF. MVMT (VEH) 

14  CALL  ELEV  (X. CURRENT  (VEH)  ,Y. CURRENT  (VEH) )  YIELDING  Z. CURRENT  (VEH) 

15  LOOP 

16  RETURN 

17  END 


Figure  6.     Routine  INIT.POS 
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B.     Routine  BEST.POS 

The  BEST.POS   routine  selects  the  best  (first)   available  position 
number  for  a  vehicle  to  occupy  in  the  platoon's  position  area. 


Input  Argument 
VEH 

Local  Variables 

ELEM 


Global  Variables  Used 
none 

Vehicle  Attributes  Used 
PLT 
POS.IN.PLT.AREA 


Routines  Called 
None 

Events  Scheduled 
None 

Sets  Used 

PLT. UNIT 

Code  -  see  Figure  7. 

Brief  Description 

Line  6 

Lines  8-13 


pointer  of  the  vehicle  to  be  positioned 


pointer  to  other  vehicles  in  VEH's 
platoon 

the  position  number 


platoon  number 

position  number  to  be  selected  by  the 
routine,  also  position  number  for 
other  vehicles 


the  platoon  to  which  the  vehicle  belongs 


Lines   14-15 


considers  J  =  1,2,...,  number  of  vehicles  in  this  platoon 

for  each  such     J,   loop  over  all   elements  in  the  platoon. 
If  any  element  has  a  POS.IN.PLT.AREA  =  J,   then  position 
J     is  already  occupied  and  thus  not  available  for  this 
vehicle.     Hence  go  to   'NEXT.POS'    to  try  the  next     J     value 

If  position     J     is  available—use  it  and  return. 
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1  ROUTINE  FOR  BEST. PCS  GIVEN  VEH 

2  "  CALLED  WHEN  VEHICLE  REACHES  END  OF  MOVEMENT  ROUTE  CLOSE  TO  NEW 

3  "  POSITION  AREA.   CHOOSES  BEST  EMPTY  POSITION  IN  HIS  PLATOON  AREA 
14  "  FOR  HIM  TO  OCCUPY. 

5  DEFINE  VEH,  J.  ELEM  AS  INTEGER  VARIABLES 

6  FOR  J  =  1  TO  N. PLT. UNIT  (PLT  (VEH) )   "NO.  ELEMENTS  IN  PLT  SET  TO  WHICH  VEH  BELONG 

7  DO 

8  FOR  EACH  ELEM  IN  PLT. UNIT  (PLT  (VEH) ) 

9  DO 

10  IF  POS. IN. PLT. AREA  (ELEM)  EQUALS  J  "POSITION  J  IS  ALREADY  FULL 

11  GO  TO  NEXT. POS 

12  ELSE 

13  LOOP   "  TO  SEE  IF  NEXT  ELEMENT  OF  PLT  IS  IN  POS  J 

14  LET  POS. IN. PLT. AREA  (VEH)  =  J 

15  RETURN 

16  'NEXT. POS' 

17  LOOP    "BACK  TO  TRY  NEXT  BEST  POSITION 

18  END 


Figure  7.  Routine  BEST. POS 
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Use 

Enter  with  vehicle  pointer. 

On  entry,  POS. IN.PLT.AREA(VEH)  should  be  zero. 

On  exit  POS. IN. PLT. AREA  will  have  been  set  to  an  integer  J 
which  is  the  first  available  position  number.  If  the  positions 
are  stored  in  prioritized  order,  the  first  available  position 
will  be  the  best  available. 

C.   Routine  RT.SEL 

The  RT.SEL  routine  selects  the  route  to  be  used  for  a  given  move- 
ment and  sets  the  ROUTE,  NEXT.MCP,  and  DIR.ON.RT  attributes  of  the  vehicle  to 
correspond  to  this  route: 


Input 

Variables 

VEH 

Local 

Variables 

Al 

A2 

P 

I,  J 

AREA 

Global 

i   Variables  Us 

ed 

MOVE. DATA 

Vehic' 

le  Attributes 

Used 

PLT 

AREA. START 

AREA. END 

ROUTE 

DIR.ON.RT 

NEXT.MCP 

Routii 

les  Called 

pointer  to  vehicle 

position  area  numbers  defining 
the  movement  desired 
platoon  number 

array  subscript 
position  area  number 


input 


output 


None 
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Events  Scheduled 
None 

Code  -  see  Figure  8. 

Brief  Description 

Lines  5-6  define  (A1,A2)   as  the  area  number  pair  with  Al   <  A2 

Lines  9-15  search  the  MOVE. DATA  array  to  find  the  pair  (A1,A2) 

Line  16  sets   ROUTE 

Lines   17-20  set  NEXT.MCP  and  DIR.ON.RT  to  indicate  forward 

movement  along  Route 

Lines  21-25  set  NEXT.MCP  and  DIR.ON.RT  to  indicate  backward 

movement  along  Route 

Lines  8  and  27       if  no  match  is  found  to  (A1,A2)  the  ROUTE  =0  on 


return. 


Use 


RT.SEL  should  only  be  used  once  at  the  start  of  a  movement 
between  areas  (MV. STATE  =  1)  because  it  sets  NEXT.MCP  to  the 
closest  end  of  the  route  regardless  of  the  actual  position 
of  the  vehicle. 

It  is  called  automatically  in  the  MOVE  routine  whenever 

MV. STATE  =  1. 

On  entry,  the  AREA. START  and  AREA. END  attributes  completely 

define  the  desired  move. 

On  exit,  the  ROUTE,  NEXT.MCP,  DIR.ON.RT  attributes  have  been 

set  to  define  details  of  the  desired  move. 

D.   ROUTINE  MOVE. LIMITS 

The  MOVE. LIMITS  subroutine  is  responsible  for  determining  limits 
on  the  speed  and  acceleration  with  which  the  vehicle  can  move.  Several  dif- 
ferent versions  of  this  routine  can  be  written,  depending  on  the  degree  of 
resolution  desired  for  movement. 
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1  ROUTINE  FOR  RT.SEL  GIVEN  VEH 

2  "  CALLED  WHEN  VEHICLE  FIRST  LEAVES  fi  POSITION  AREA.   CHOOSES  THE  ROUTE 

3  "  ALONG  WHICH  VEH  WILL  MOVE  TO  REACH  NEXT  POSITION  AREA 
U  DEFINE  VEH,A1.A2,P,I.J.AREA  AS  INTEGER  VARIABLES 

5  LET  Al  =  MIN.F (AREA. START (VEH) .AREA. END  (VEH)) 

6  LET  A2  =  MAX. F  (AREA. START  (VEH)  .AREA. END  (VEH)) 

7  LET  P  =  PLT  (VEH) 

8  LET  ROUTE  (VEH)  =  0 

9  FOR  1  =  1  TO  D1M.F  (MOVE. DATA  (P. *) ) /3 

10  DO 

11  LET  J  =  3x1 

12  LET  AREA  =  MOVE. DATA  (P. J-2) 

13  IF  AREA  IS  GREATER  THAN  Al  RETURN   "NO  MATCH  FOR  AREA  NUMBERS  FOUND 

14  ELSE 

15  IF  AREA  EQUALS  Al  AND  MOVE. DATA  (P, J-l )  EQUALS  A2 

16  LET  ROUTE  (VEH)  =  MOVE. DATA  (P. J) 

17  IF  AREA. START  (VEH)  IS  LESS  THAN  AREA. END  (VEH)   "NORMAL  DIRECTION 

18  LET  DIR.ON.RT  (VEH)  =  0 

19  LET  NEXT. MCP  (VEH)  =  1 

20  RETURN      "WITH  FORWARD  ROUTE 

21  ELSE 

22  LET  DIR.ON.RT  (VEH)  =  1 

23  LET  NEXT. MCP  (VEH)  =  DIM. F  (ROUTE. DATA  (ROUTE  (VEH) , *) ) /3 
21  RETURN      "WITH  BACKWARD  ROUTE 

25  ELSE 

26  LOOP      "BACK  TO  TRY  NEXT  SEGMENT  OF  MOVE. DATA  ARRAY 

27  RETURN    "WITH  NO  ROUTE  --  MATCH  NOT  FOUND 

28  END 


Figure  8.     Routine  RT.SEL 
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A  high  resolution  movement  model  might  consider  digitized  limiting 
speed  and  acceleration  values  based  on  extensive  terrain  analysis  and  on 
detailed  vehicle  characteristics. 

A  low  resolution  movement  model  might  set  these  values  constant  for 
each  vehicle  type  independent  of  terrain  details;  other  resolutions  are  not 
difficult  to  imagine. 

In  any  case,  the  influence  of  terrain  on  movement  can  be  concentrated 
in  the  MOVE. LIMITS  routine  thus  enhancing  the  flexibility  of  the  model.  In 
this  report  we  will  not  detail  the  MOVE. LIMITS  routine  as  it  depends  so  much 
on  the  terrain  representation  chosen.  The  general  characteri sites  of  the 
routine  are,  however, 
Input  Variables 

Pointer  to  the  vehicle 

Dimensionless  terrain  slope(rise  *  run) 


VEH 
SLOPE 
Output  Variables 
SPEED 


ACCEL 


The  target  speed  which  the  movement  model  should  try 
to  match  in  this  move  increment.  This  speed  may  de- 
pend on  terrain,  the  unit's  desired  maneuver  speed, 
obstacles  or  minefields,  etc.  In  particular,  if  the 
vehicle's  MV. STATE  is  3,  SPEED  should  be  set  to  0  (Stop) 

(>  0)  limit  on  allowed  acceleration  of  the  vehicle. 
Again  this  may  be  modelled  in  varying  degrees  of  detail. 

(<  0)  limit  on  allowed  deceleration. 


DECEL 
Calling  Sequence 

CALL  MOVE.LIMITS(VEH,  SLOPE)  YIELDING  SPEED,  ACCEL,  DECEL 
Called  From 

MOVE  at  start  of  each  movement  increment 
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E.   Routine  ELEVG 

Routine  ELEVG  provides  the  macro  terrain  representation  for  the 
simulation.  For  any  X,  Y  coordinates  on  the  battlefield  it  provides  the 
Z  (vertical  elevation)  coordinate  and  the  gradient  components  GX  and  GY 
of  the  terrain.  ELEVG  is  similar  to  MOVE. LIMITS  in  that  it  is  independent 
of  the  rest  of  the  move  model  and  can  be  realized  in  several  different  ways 
depending  on  the  study  requirements.  For  example,  tabletop  terrain  is  obtained 
by  having  ELEVG  return  a  constant  elevation  and  zero  gradient.  The  current 
STAR  implementation  utilizes  functionally   coded  terrain,  and  an  ELEVG  routine 
has  been  written  for  that  representation.  Changing  to  a  digitized  terrain 
representation  merely  requires  that  the  appropriate  routine  be  named  ELEVG 
and  dropped  into  place.  It  is  important  to  note  that  this  (and  any  other) 
routine  may  be  written  in  either  FORTRAN  or  SIMSCRIPT. 


Input 

Variables 

X 

1 

Y 

Output  Variables 

Z 

GX 

GY 

} 

} 


map  coordinates 

elevation 

gradient  components 


Attributes  Used 
None 

Routines  Called 
None 

Events  Scheduled 
None 

T21 
Code   See  other  STAR  publications    for  the  currently  used  functional 

terrain  model  and  the  corresponding  ELEVG  program. 
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F.  Routine  MOVE 

The  MOVE  routine  updates  the  location,  direction,  and  speed  of  a 
vehicle  to  the  current  simulation  time.  The  routine  is  rather  long,  but  not 
overwhelmingly  complex.  It  is  useful  to  consider  it  as  five  sequential  parts 
which  are  almost  always  executed  in  strict  sequence: 
i)   compute  a  destination  point 
ii)   compute  direction  and  angles  from  current  location 

to  destination  point 
iii)   relate  distance,  time,  speed,  and  acceleration  to 

define  the  move  increment 
iv)   update  location  and  time  for  this  move  increment 
v)   check  whether  finished.  If  not,  go  back  to  i)  or  iii). 
Parts  i)  and  iii)  comprise  the  bulk  of  the  code  because  each  contains 
several  alternative  possible  computational  sequences.  Part  i)  determines 
the  various  modes  of  movement  as  listed  earlier  in  this  report,  and  Part  iii) 
must  consider  move  increments  limited  by  either  time  or  distance  for  any  given 
speed  and  acceleration  limits.  In  describing  the  logic  of  the  move  routine  we 
will  use  this  breakdown  by  the  5  listed  program  segments. 


Input  Variable 
VEH 

Local  Variables 
numerous 


Global  Variables  and  Arrays 
POSITION 
ROUTE. DATA 


all  are  listed  in  define  statements 
at  beginning  of  program--we  explain 
them  as  appropriate  in  discussing 
program  logic. 

positions  in  areas 
MCP  coordinates 
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FORM. OFFSET  formations 

FOR.CHG.INT  distance  within  which  formation 

changes  are  accomplished 

MAX.DIST. INCR  max  distance  to  move  before  re- 

evaluating terrain. 

Routines  Called 

RT.SEL 
BEST.POS 
ELEVG 
MOVE. LIMITS 

Events  Scheduled 
None 

Code 

We  break  the  code  into  several   segments  and  discuss  each  in  turn. 

Segment  0)     Declarations,   Initialization,  and  MV. STATE  filters  -  see  Fig.    9. 

Lines  19-22  if  MV. STATE  is  1  we  do  a  route  select  (and  set 

MV. STATE  =  2  for  all   future  calls  to  move  along  this   route) 

Lines  23         the  time  interval   over  which  this  vehicle's  movement  is 

to  be  computed  runs  from  T.SPD  to  TIME.V  (current  simulation  time) 

Segment  i)     Compute  a  Destination  Point 

The  various  modes  of  movement  are  computed  in  this  segment  of  the 
code.     The  result  of  the  computation  is  a  destination  point     X.DEST,  Y.DEST 
toward  which  the  vehicle  will   move. 

The  flowchart  of  Figure  10  describes  the  possible  alternatives. 
The  Code,  given  in  Figure  12,  is  quite  long  because  of  the  several  distinct 
alternatives. 
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1  ROUTINE  TO  MOVE  GIVEN  VEH 

2  DEFINE  K  AS  AN  INTEGER  VARIABLE 

3  DEFINE  SLOPE  AS  A  REAL  VARIABLE 

4  DEFINE  REM. MOVE. TIME.  DEL.X.  DEL.Y.  D.TO.MCP.  ALPH.   SALPH. 

5  CALPH.  GRAD.X.  GRAD.T.  SPD. LIMIT.  ACCEL. LIMIT.  DECEL. LIMIT, 

6  DIST. LIMIT.  DEL. SPD.  DIST.INCR.  T1ME.INCR  AS  REAL  VARIABLES 

7  DEFINE  DIST.REQ.  TIME. RED,  ZERO. LEVEL  AS  REAL  VARIABLES 

8  DEFINE  X.DEST,  T.DEST,  DIR.  CX,  CT.  NX.  NY.  LX.  LY.  NLX.  NLY.  PX.PY, 

9  NPX.  NPY.  X.OFF.  Y.OFF.  D.TO.DEST  AS  REAL  VARIABLES 

10  DEFINE  THETA,  CTH,  STH  AS  REAL  VARIABLES 

11  DEFINE  VEH.  FINAL  AS  INTEGER  VARIABLES 

12  DEFINE  MST.  RT.  NM,  MCP.INC.  LM,  MCP.  D.ON.RT  AS  INTEGER  VARIABLES 

13  DEFINE  FAKE. MCP  AS  AN  INTEGER  VARIABLE 
m  DEFINE  I,  J  AS  INTEGER  VARIABLES 

15  LET  ZERO. LEVEL  =  1.0     LET  FINAL  =  0 

16  LET  MST  =  MV. STATE  (VEH) 

17  IF  MST  EQ  0  OR  MST  GE  4  RETURN 

18  ELSE 

19  IF  MST  EQUALS  1 

20  CALL  RT.SEL  (VEH) 

21  LET  MV. STATE  (VEH)  =  2 

22  ALWAYS 

23  LET  REM. MOVE. TIME  =  TIME.V  -  T. SPD  (VEH) 


Figure  9.  MOVE  Routine,  Segment  0 
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Go  to  ANGLES  (in  segment  ii) 
to  continue  moving  in  previous 
DIR.OF.MVMT  (line  27) 


User  has  swapped  AREAS  to  call  for 
reversal.  Change  DIR.ON.RT,  NEXT.MCP 
to  reflect  this  reversal  (lines  29-52) 


Have  reached  end  of  route  -- 
move  to  a  position  in  AREA. END 
(lines  54-74) 


NO 


NO 


Choose  designated 
position  as 
destination 


Choose  best 
position  as 
destination 


Figure  10. 

Flowchart  of  Segment  1 
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HAVE  WE 
ALREADY 
ARRIVED 
sJHERE 

Set  MV. STATE 
S.YES  w&  terminate 
'        r  (lines  68-73) 

N.    ?     > 

NO 
GO 

TO 

DIRN.COMP 

[in 

segment  ii) 

NO 


no  formation  offset— let 
t  destination  be  the  next 

MCP  on  the  route  (lines  78-82). 


1 


go    to  DIRN.COMP 
in  segment  ii) 


compute  destination  based  on 
first  and  second  MCP's 
(lines    85-89) 


compute  destination  based  on 
last  and  next  to  last  MCP's 
(lines  93-97) 


setup  destination  between  the  two 
MCP's  to  get  onto  route  and  into 
formation  (lines  101-114) 


I 


go  to  DIRN.COMP 
(in  segment  ii) 


Figure.   10  (Continued) 
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INTERMED' 


Setup  last  &  next  MCP's 
for 

moving  along 
intermediate 
portion  of  route 
(lines  116-128) 


YES 


use  formation  code 
■^         on  the  vehicle 
(line   130) 


use  formation  code 
on  the  route 
(line   129) 


let  destination  be  a  point  FOR.CHG.INT  ahead  along  the  route  and 
offset  by  the  desired  formation  offsets   (lines   132-144) 

if  destination  is  closer  than  next  MCP,  call   it  a  fake  MCP  to  avoid 
messing  up  MCP  count  later  in  the  program  (lines   145-148 


go  to  DIRN.COMP 
(in  segment  ii ) 


Figure  10.  (Continued) 


The  destination  point  computations  are  a  bit  intricate,  especially 
in  lines  116-148   Reference  to  Figure  11  will  clarify  the  several  intermed- 
iate variables  used. 

In  each  case,except  where  RT=  05  we  leave  this  code  segment  with  a 
destination  point  X.DEST,  Y.DEST  and  a  distance  D. TO. MCP  which  is  the  distance 
we  can  move  before  having  to  recompute  a  new  destination  point.  The 
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NEXT.MCP 
(NX, NY) 


FOR.CHG.INT 


computed 

destination 

point 

(X.DEST,   Y.DEST) 

X.OFF 


current 
location  of 
vehicle^ 
(CX,CY) 


(PX,PY)  projection  of  current 

location  onto  route  segment 


(LX.LY)  previous  MCP 


Route 


Figure  11.  Geometry  of  Destination  Computation 

in  Intermediate  Case 
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24  LET  RT  =  ROUTE  (VEH) 

25  LET  D.ON.RT  =  DIR. ON. RT  (VEH) 

26  LET  FAKE. MCP  =  0 

27  IF  RT  EQURLS  0    LET  O.TO.MCP  =  RINF.C    GO  TO  ANGLES 

28  ELSE 

29  "CONSISTENCY  CHECK  FOR  POSSIBLE  TURNAROUND 

30  IF  AREA. START (VEH)  LT  AREA. END  (VEH) 

31  IF  D.ON.RT  EQ  0   GO  TO  NEW.MCP 

32  ELSE  LET  D.ON.RT  =  0 

33  LET  DIR.ON.RT(VEH)  =  0 

34  LET  K  =  DIM.F  (ROUTE. DATA (RT.x)  ) /3 

35  LET  MCP  =  NEXT. MCP (VEH) 

36  IF  MCP  EQ  0  LET  NEXT. MCP  (VEH)  =  1 

37  ELSE  IF  MCP  EQ  K   LET  NEXT. MCP  (VEH)  =  0 

38  ELSE  LET  NEXT.  MCP (VEH)  =  MCP  +  1 

39  ALWAYS 

40  ALWAYS 

41  ELSE  "AREA. START  GT  AREA. END 

42  IF  D.ON.RT  EQ  1  GO  TO  NEW.MCP 

43  ELSE  LET  D.ON.RT  =  1 

44  LET  DIR.ON.RT  (VEH)  =  1 

45  LET  K  =  DIM. F  (ROUTE. DATA (RT.*)) /3 

46  LET  MCP  =  NEXT. MCP (VEH) 

47  IF  MCP  EQ  0  LET  NEXT. MCP  (VEH)  =  K 

48  ELSE  IF  MCP  EQ  1  LET  NEXT. MCP  (VEH)  =  0 

49  ELSE  LET  NEXT. MCP (VEH)  =  MCP  -1 

50  ALWAYS 

51  ALWAYS 

52  ALWAYS 

53  'NEW.MCP*  LET  MCP  =  NEXT. MCP (VEH)  LET  NM  =  MCP  *  3 

54  IF  MCP  EQUALS  0   "MOVE  TO  POSITION  IN  AREA. END 

55  IF  POS.IN.PLT.AREA(VEH)  EQUALS  0,  CALL  BEST. POS (VEH)   "SETTING  POS.IN.PLT.A 

56  ALWAYS 

57  LET  1  =  PLT(VEH)  LET  K  =  POS. IN. PLT. AREA  (VEH)  *  3 

58  FOR  J  =  1  TO  DIM.F(POSITIQN(J.m.*))  WITH  POSITION  (I . J. 1)  EQUALS 

59  AREA. END  (VEH) 

60  DO 

61  LET  X.DEST  =  POSITION  (I , J, K-l) 

62  LET  Y.DEST  =  POSITION (I , J, K) 

63  LET  DIR  =  POSITION  (I . J. K+l) 

64  LOOP 

65  LET  D. TO. MCP  =  SQRT. F  (  (X. DEST-X. CURRENT  (VEH) ) **2  + 

66  (Y. DEST-Y. CURRENT  (VEH) ) **2) 

67  IF  D. TO. MCP  LESS  THAN  ZERO. LEVEL, 

68  LET  MV. STATE  (VEH)  =  4 

69  LET  DIR.OF.MVMT  (VEH)  =  DIR 

70  LET  PRI.DIR(VEH)  =  DIR 

71  LET  SPD(VEH)  =  0. 

72  LET  FINAL  =  1 

73  GO  TO  NEW.INCR 

Figure  12.  MOVE  Routine,  Segment  i 
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74  ELSE 

75  60  TO  DIRN.COMP 

76  ELSE 

77  IF  FORM. CODE  (VEH)  EQUALS  0       "GO  DIRECTLY  TO  NEXT  MCP 

78  LET  X.DEST  =  ROUTE. DATA  (RT.NM-2) 

79  LET  T.DEST  =  ROUTE. DATA  (RT, NM-1) 

80  LET  D. TO. MCP  =  SORT. F  (  (X.DEST-X. CURRENT (VEH) )  **2  ♦ 

81  (Y. DEST-T. CURRENT (VEH)  )  **2) 

82  GO  TO  DIRN.COMP 

83  ELSE      "MOVE  ALONG  ROUTE  OFFSET  BY  FORMATION 

84  IF  MCP  EQUALS  1  AND  D.ON.RT  EQUALS  0   "TOWARD  FIRST  MCP 

85  LET  NX  =  ROUTE. DATA  (RT.4) 

86  LET  NY  =  ROUTE. DATA  (RT. 5) 

87  LET  LX  =  ROUTE. DATA  (RT, 1) 

88  LET  LY  =  ROUTE. DATA  (RT, 2) 

89  LET  I  =  ROUTE. DATA  (RT, 3) 

90  ELSE 

91  LET  K  =  OIM.F(ROUTE.DATA(RT.k)) 

92  IF  MCP  EQUALS  K/3  AND  D.ON.RT  EQUALS  1    "TOWARD  LAST  MCP  GOING  BACKWARD 

93  LET  NX  =  ROUTE. DATA  (RT.K-5) 

94  LET  NY  =  ROUTE.  DATA (RT,  K-4) 

95  LET  LX  =  ROUTE. DATA  (RT.K-2) 

96  LET  LY  =  ROUTE. DATA  (RT.K-1) 

97  LET  I  =  ROUTE. DATA  (RT.K-3) 

98  ELSE  GO  TO  INTERMED 

99  ALWAYS 

100  ALWAYS 

101  LET  NLX  =»  NX-LX       LET  NLY  =  NY-LY 

102  IF  I  EQUALS  0,   LET  I  =  FORM. CODE  (VEH) 

103  ALWAYS 

104  LET  J  *     FORM. POS  (VEH)  *  2 

105  LET  X.OFF  =  FORM. OFFSET  (1 , J-l) 

106  LET  Y.OFF  =  FORM. OFFSET (I . J) 

107  LET  THETA  =  ARCTAN. F  (NLY, NLX) 

108  LET  CTH  =  COS. F  (THETA) 

109  LET  STH  «  SIN.F  (THETA) 

110  LET  X.DEST  *  LX  +  (X.OFF  ♦  FOR.CHG. INT) *CTH  -  Y.0FF*STH 

111  LET  Y.DEST  =  LY  ♦  (X.OFF  ♦  FOR. CHG. INT) *STH  +  Y.OFF  *  CTH 

112  LET  Q. TO. MCP  =  SORT. F ( (X.DEST-X. CURRENT (VEH) ) «*2   + (Y.DEST-Y. CURRENT (VEH) ) 

113  km2) 

114  GO  TO  DIRN.COMP 

115  'INTERMED'      "TO  HERE  FOR  INTERMEDIATE  MCP'S  ON  ROUTE 

116  LET  CX  =»  X. CURRENT  (VEH)      LET  CY  -  Y. CURRENT  (VEH) 

117  IF  D.ON.RT  EQUALS  0   LET  LM  =  NM  -  3 

118  ELSE  LET  LM  =  NM  ♦  3 

119  ALWAYS 

120  LET  NX  =  ROUTE. DATA (RT.NM-2) 

121  LET  NT  =  ROUTE. DATA (RT. NM-1) 

122  LET  LX  =  ROUTE. DATA (RT.LM-2) 

123  LET  LY  =  ROUTE. DATA  (RT.LM-1) 


Figure  12.   (Continued) 
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124  LET  NLX  =  NX  -  LX  LET  NLT  =  NY  -  LY 

125  LET  flLPH  =- ( (CX-NX) *NLX  +  (CY-NY)  *NLY)  /  (NLXmNLX  +  NLY*NLY) 

126  LET  PX  =  flLPH  *  LX  +  tl.  -  flLPHJ  *  NX 

127  LET  PT  =  flLPH  *    LY  +  (1.  -  flLPH)  *  NY 

128  LET  NPX  =  NX  -  PX  LET  NPY  =  NY  -  PY 

129  LET  I  =  ROUTE.  DATA  (RT.NM  +  3* (D.  ON.  RT-1)) 

130  IF  I  EQUALS  0.   LET  I  =  FORM. CODE  (VEH) 

131  ALWAYS 

132  LET  J  -     FORM. POS (VEH)  *    2 

133  LET  X.OFF  =  FORM. OFFSET (1 . J-l) 

134  LET  Y.OFF  =  FORM. OFFSET (I , J) 

135  LET  D.TO.MCP  =  SORT. F  (NPX*NPX  +  NPY*NPY) 

136  IF  D.TO.MCP  LESS  THAN  ZERO. LEVEL    GO  TO  MCP. REACHED 

137  ELSE 

138  LET  THETA  =  ARCTAN. F  (NLY, NLX) 

139  LET  CTH  =  COS. F  (THETA) 
110  LET  STH  =  SIN. F  (THETA) 

141  LET  ALPH  =  FQR.CHG.INT  /  D.TO.MCP 

142  LET  X.DEST  =  PX  +  flLPH  *    NPX  ♦  X.OFF   *  CTH  -  Y.OFF  x  STH 

143  LET  Y.DEST  =  PY  +  flLPH  «  NPY  ♦  Y.OFF  *  CTH  ♦  X.OFF  *  STH 

144  LET  D.TO.DEST  =  SORT. F  (  (X. DEST-CX) **2  +  (Y.DEST  -CY) **2) 

145  IF  D.TO.DEST  IS  LESS  THAN  D.TO.MCP 

146  LET  D.TO.MCP  =»  D.TO.DEST 

147  LET  FAKE. MCP  =  1 

148  ELSE  LET  FAKE. MCP  =  0 

149  ALWAYS 


Figure  12.   (Continued) 
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destination  point  may  be  an  MCP,  a  position  in  a  position  area,  or,   in 
the  case  of  movement  in  formation,  a  so-called  fake-MCP,  created  tem- 
porarily to  get  the  vehicle  into  the  proper  formation  position. 
Segment  ii)  Compute  direction  and  angles 

The  brief  code  segment  given  in  Figure  13  determines  the 
DIR.OF.MVMT  to  the  destination  point  and  computes  some  trig  functions  for 
later  use  in  segment  iv). 
Segment  iii )         Relate  distance,  time,  speed  and  acceleration. 

Segment  iii   compares  the  remaining  move  time  to  the  distance 
to  be  moved  in  this  increment.     Depending  on  the  desired  target  speed, 
acceleration  capabilities,  and  the  time  and  distance  limits,  a  DIST.INCR,  a 
TIME. INCR,  and  a  final   SPD  are  computed   (see  Figure  14). 
Lines  159-165       Sample  the  terrain  using  routines  ELEVG  and 

MOVE.   LIMITS  to  get  a  target  speed  SPD. LIMIT, 

and  limits  ACCEL. LIMIT  >  0  and  DECEL. LIMIT  <  0. 
Line  166  Computes  a  distance  limit  for  the  move.     Note  the 

user  supplied  MAX. DIST.INCR  which  forces  periodic 

re-sampling  of  the  terrain. 
Lines  169-178   Handle  the  frequently  occurring  special  case  where 

speed  is  constant  throughout  the  move  increment. 

Simple  manipulation  of  the  movement  equation 

d  ■  V 

yields  a  time  increment  and  a  distance  increment 
for  the  move. 
Lines  179-205   Consider  the  more  complex  situation  where  acceleration 

occurs.  The  movement  equation 

1    2 
d  =  vQt  +  j  at  c 
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150  'DIRN.COMP' 

151  IF  D.TO.MCP  15  LESS  THAN  ZERO. LEVEL, 

152  GO  TO  MCP. REACHED 

153  ELSE 

154  LET  DEL.X  *  X.DEST  -  X. CURRENT  (VEH) 

155  LET  DEL.Y  =  Y.DEST  -  Y. CURRENT  (VEH) 

156  LET  DIR. OF. MVMT  (VEH)  =  ARCTAN. F  (DEL. Y, DEL. X) 

157  'ANGLES' 

158  LET  SALPH  =  SIN. F (DIR. OF. MVMT  (VEH) )      LET  CALPH  =  COS. F (DIR. OF. MVMT (VEH)  ) 


Figure  13.     MOVE  Routine  -  Segment  ii 
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159  'NEW. INCR'    CALL  ELEVG  GIVEN  X. CURRENT  (VEH) ,  Y. CURRENT  (VEH)  YIELDING 

160  Z. CURRENT  (VEH) ,  GRAD.X.  GRAO.Y 

161  IF  FINAL  EQUALS  1,  GO  TO  END. MOVE 

162  ELSE 

163  LET  SLOPE  =  GRAD.X  «  CflLPH  ♦  GRAD.Y  x  SALPH 

161  CALL  MOVE. LIMITS  GIVEN  VEH.  SLOPE          YIELDING  SPD. LIMIT.  ACCEL. LIMIT. 

165  DECEL. LIMIT 

166  LET  DIST. LIMIT  =  MIN. F  (0. TO. MCP.  MAX.DIST. INCR) 

167  LET  DEL.SPD=  SPD. LIMIT  -  SPD  (VEH) 

168  IF  ABS.F(DEL.SPD)  IS  LESS  THAN  0.1 

169  "EAST  CASE  --  NO  ACCELERATION  -- 

170  LET  DIST. INCR  =  REM. MOVE. TIME  *   SPD. LIMIT 

171  IF  DIST. INCR  IS  GREATER  THAN  DIST. LIMIT. 

172  "MOVE  STOPPED  BY  DISTANCE  LIMIT 

173  LET  DIST. INCR  =  DIST. LIMIT 

174  LET  TIME. INCR  =  DIST. INCR  /  SPD. LIMIT 

175  ELSE 

176  "MOVE  STOPPED  BY  TIME  LIMIT 

177  LET  TIME. INCR  =  REM. MOVE. TIME 

178  ALWAYS  GO  TO  MOVE. IT 

179  ELSE   "HARD  CASE  --  ACCELERATION  REQUIRED  — 

180  IF  DEL. SPD  IS  LESS  THAN  0.  LET  ACCEL. LIMIT  =  DECEL. LIMIT 

181  ALWAYS  LET  TIME.REQ  =  DEL. SPD  /  ACCEL. LIMIT 

182  LET  DIST.REQ  =  SPO (VEH) *TIME.REQ  +  0.5  *  ACCEL. LIMIT  *  TIME.REQ  **2 

183  IF  TIME.REQ  IS  GREATER  THAN  REM. MOVE. TIME. 

184  "SPD. LIMIT  CANNOTBE  ATTAINED  IN  REMAINING  TIME,  SO  CHANGE  LIMIT 

185  LET  SPD. LIMIT  *  SPD  (VEH)  ♦  ACCEL. LIMIT  x  REM. MOVE. TIME 

186  LET  DIST. INCR  *  SPD  (VEH)  *  REM. MOVE. TIME  ♦  0.5  *  ACCEL. LIMIT  * 

187  REM. MOVE. TIME  **  2 

188  ELSE    "SPO. LIMIT  CAN  BE  ATTAINED 

189  LET  DIST. INCR  *  DIST.REQ  ♦  (REM. MOVE. TIME  -  TIME.REQ) «SP0. LIMIT 

190  ALWAYS 

191  IF  DIST. INCR  IS  LESS  THAN  DIST. LIMIT. 

192  "MOVE  WILL  BE  STOPPED  BY  TIME  LIMIT 

193  LET  TIME. INCR  =  REM. MOVE. TIME 

194  LET  SPD  (VEH)  =•  SPD. LIMIT 

195  ELSE     "MOVE  STOPPEO  BY  DISTANCE  LIMIT 

196  LET  DIST. INCR  =  DIST. LIMIT 

197  IF  DIST. LIMIT  IS  LESS  THAN  DIST.REQ. 

198  LET  TIME. INCR  =  ISQRT. F  (SPD  (VEH) **2*2.xACCEL.LlMIT*0IST. LIMIT) 

199  -SPD  (VEH)) /ACCEL. LIMIT 

200  ADD  TIME. INCR  *  ACCEL. LIMIT  TO  SPD  (VEH) 

201  ELSE 

202  LET  TIME. INCR  =  TIME.REQ  ♦ (DIST. LIMIT-DIST.  REQ)  /SPD. LIMIT 

203  LET  SPD  (VEH) =SPD. LIMIT 

204  ALWAYS 

205  ALWAYS 


Figure  14.  MOVE  Routine,  Segment  iii 
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is  manipulated  in  various  ways  to  again  yield  a  time 
and  a  distance  increment  and  also  the  final  speed  of 
the  vehicle  at  the  end  of  the  move  increment. 

In  each  case,  the  results  of  segment  iii),  TIME.INCR,  DIST.INCR,  and  SPD(VEH) 

are  passed  on  to  segment  iv)  to  actually  perform  the  movement  increment. 

Segment  iv)     Update  location  and  time 

Segment  tv)  performs  the  move  by  changing  the  vehicle's 

X  and  Y  coordinates.  The  code  which  follows  is  self  explanatory. 

206  'MOVE. IT'      SUBTRACT  TIME.INCR  FROM  REM. MOVE. T IME 

207  ADD  DIST.INCR  *  CALPH  TO  X. CURRENT  (VEH) 

208  ADD  DIST.INCR  *  SALPH  TO  Y. CURRENT  (VEH) 

209  SUBTRACT  DIST.INCR  FROM  D.TO.MCP 

Segment  v)  Check  whether  finished  (see  Figure  15  for  Code). 

Various  occurrences  can  end  a  MOVE  call.   If  the  move  time  has 
expired,  then  except  for  updating  the  elevation  we  are  finished  (line  210). 

If  D.TO.MCP  has  been  exceeded,  then  a  new  direction  computation  is 
needed,  and,  if  the  MCP  reached  is  a  real  one,  we  want  to  aim  toward  the  next 
MCP  (lines  213-233). 

If  a  mine  detonation  or  a  minefield  entry  or  exit  has  occurred,  we 
must  pause  to  assess  the  implications  (damage,  lower  mine  plow,  ...).  Coding 
for  these  functions  is  not  yet  written,  but  the  tests  for  them  will  be  in- 
cluded at  the  beginning  of  this  segment. 

At  the  end  of  the  move  time  specified  for  this  move,  the  T.SPD(VEH) 
is  updated  to  the  current  simulation  time  and  control  returns  to  the  calling 
program,   (lines  235-236). 
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210 

IF  REM. MOVE. TIME  IS  LESS  THAN   0.01      LET  FINAL  =  1 

211 

REGARDLESS 

212 

IF  D.TO.MCP  IS  LESS  THAN  ZERO. LEVEL, 

213 

'MCP. REACHED' 

2m 

IF  FAKE. MCP  EQUALS  1 

215 

LET  FAKE. MCP  =  0 

216 

GO  TO  NEW. MCP 

217 

ELSE 

218 

IF  MCP  =  0 

219 

LET  FINAL  *  1 

220 

GO  TO  NEW. MCP 

221 

ELSE 

222 

IF  DIR.ON.RT  (VEH)  EQUALS  0    "MCP  NUMBERS  INCREASE 

223 

IF  NEXT.MCP(VEH)  EQUALS  DIM. F  (ROUTE. DATA  (ROUTE  (VEH 

224 

LET  NEXT. MCP  (VEH)  =  0 

225 

GO  TO  NEW. MCP 

226 

ELSE 

227 

ADD  1  TO  NEXT. MCP  (VEH)      GO  TO  NEW. MCP 

228 

ELSE      "MCP  NUMBERS  DECREASE 

229 

IF  NEXT. MCP  (VEH)  EQUALS  1 

230 

LET  NEXT. MCP (VEH) =0 

231 

GO  TO  NEW. MCP 

232 

ELSE 

233 

SUBTRACT  1  FROM  NEXT. MCP  (VEH)     GO  TO  NEW.MC 

231 

ELSE   GO  TO  NEW.INCR 

235 

'END. MOVE'    LET  T.SPD(VEH)  =  TIME.V 

236 

RETURN 

237 

END 

.*)  )  /3 


Figure  15.  MOVE  Routine,  Segment  v, 


-  44  - 


IV.   INTERFACE 

This  section  of  the  report  will  concentrate  on  the  data  flow  between 
the  MOVE  routine  and  the  rest  of  the  STAR  model,  showing  how  to  use  the  MOVE 
routine  in  each  of  the  possible  modes. 
A.  MODES  OF  USE 

1.  The  most  frequently  used  mode  will  be  movement  from  one  position 
area  to  another.  To  initiate  this  (combination)  mode  the  STAR  model  should 
set  AREA. START  and  AREA. END  to  the  appropriate  area  numbers,  set  MV. STATE  =  1, 
set  T.SPD  to  the  time  at  which  the  move  was  to  have  started,  and  call  MOVE(VEH). 
The  MOVE  routine  will  select  the  route,  start  the  vehicle  moving,  and  return 
with  a  MVSTATE  of  2.  Subsequent  calls  to  MOVE  to  continue  the  movement 

should  leave  MV.STATE=  2.  When  the  vehicle  reaches  its  final  position  in 
AREA. END,  the  MOVE  model  will  return  MV. STATE  =  4. 

If  the  user  wants  the  vehicle  to  select  the  best  position  in 
AREA. END,  POS. IN.PLT.AREA(VEH)  should  be  set  to  0  prior  to  the  first  call 
to  MOVE  and  left  unchanged  thereafter. 

Movement  will  be  in  the  formation  specified  by  the  vehicle  (or 
forced  by  the  terrain)  unless  FORM. CODE  =  0  in  which  case  movement  will  be 
along  the  route. 

Prior  to  some  subsequent  call  to  move;  (while  MV. STATE  still  =  2), 
if  AREA. START  and  AREA. END  are  interchanged,  then  the  MOVE  model  will  reverse 
the  vehicle's  direction  of  movement  along  the  route. 

2.  A  simple  mode  of  movement  is  straight  line  motion  from  current 
location  to  a  position  area  (where  current  location  may  also  be  in  a  position 
area) . 
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To  achieve  this,  the  user  should  set 

MV.STATE(VEH)  =  2  (no  route  select) 

ROUTE(VEH)  =  any  number  except  0 

(the  route  will  not  actually  be  used) 

NEXT.MCP(VEH)  =  0 

AREA.END(VEH)  =  the  desired  area 

T.SPD  (VEH)   =  the  time  at  which  the  move  was  to  have  started 
POS.IN.PLT.AREA(VEH)  =  the  desired  position  #  (or  0  for  best 

position)  and  call  MOVE(VEH). 
The  straight  line  movement  will  continue  with  subsequent  MOVE  calls 

until  the  vehicle  reaches  the  position  (signalled  by  MV. STATE  =  4  on  return). 

Formations  cannot  be  used  in  this  mode  as  there  are  no  routes. 

3.  The  simplest  mode  of  movement  is  movement  in  a  specified  direction. 
To  achieve  this  the  user  should  set 

DIR.OF.MVMT(VEH)  =  the  desired  direction 
MV. STATE (VEH)  =  2 
ROUTE(VEH)  =  0 
and  call  MOVE(VEH).  Movement  in  the  specified  direction  will  continue,  upon 
subsequent  calls  to  move,  with  no  consideration  given  to  routes,  position  areas, 
or  formations.  CAUTION:  this  mode  may  lead  to  vehicles  driving  off  the  map 
with  unpredictable  (unusually  disastrous)  results  in  the  simulation.  The  user 
can  stop  the  movement  by  setting  MV. STATE  =  0,  3,  4,  or  5  at  any  time. 

4.  Stopping.  In  any  of  the  above  modes  of  movement,  if  MOVE  is  called 
with  MV. STATE  =  3,  then  the  MOVE. LIMITS  routine  should  set  SPD. LIMIT  =  0  and 
the  vehicle  will  try  to  stop.  If  the  deceleration  rate  allows,  the  vehicle 
will  stop  with  SPD(VEH)  =  0  on  return.  If  the  movement  time  is  too  short  to 
allow  stopping,  the  vehicle  will  decelerate  as  much  as  possible,  and  upon  sub- 
sequent MOVE  calls,  will  stop.  To  start  again,  set  MV. STATE  =  2  and  call  MOVE. 
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B.  Current  Use  in  STAR 

In  the  current  version  of  the  STAR  model,  the  decisions  of  when 
and  where  to  initiate  movement  are  handled  by  the  movement  decision  logic 
and  movement  coordination  logic  which  are  documented  in  reference  [4]. 
Mode  1)  above  is  used  exclusively.  Whenever  the  model  needs  to  know  the 
location  of  a  vehicle,  the  MOVE  routine  is  then  called  to  update  its  po- 
sition. This  is  done  periodically  for  all  vehicles  in  event  STEP. TIME  and 
also  at  other  times  for  individual  vehicles  involved  in  a  firing  event  as 
either  shooter  or  target. 
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V.    EXTENSIONS 

The  primary  aspects  of  movement  in  the  STAR  model  which  are  not 
covered  in  this  report  are 

a)  deciding  when  and  where  to  move  (see  ref.  [4].) 

b)  movement  of  aircraft  (see  ref.  [3].) 

c)  interaction  of  moving  vehicles  with  minefields  and  obstacles 
of  various  kinds. 

A  generil  field  structure  has  been  developed  to  include  minefields,  ob- 
stacles, and  other  battlefield  features.  This  module  is  currently 
undergoing  tests  and  will  be  the  subject  of  a  future  report. 
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